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We characterize the different morphological phases that occur in a simple one-dimensional model 
of propagation of innovations among economic agents [X. Guardiola, et. ai, Phys. Rev E 66, 026121 
(2002)]. We show that the model can be regarded as a nonequilibrium surface growth model. This 
allows us to demonstrate the presence of a continuous roughening transition between a flat (system 
size independent fluctuations) and a rough phase (system size dependent fluctuations). Finite- 
size scaling studies at the transition strongly suggest that the dynamic critical transition does not 
belong to directed percolation and, in fact, critical exponents do not seem to fit in any of the known 
universality classes of nonequilibrium phase transitions. Finally, we present an explanation for the 
occurrence of the roughening transition and argue that avalanche driven dynamics is responsible for 
the novel critical behavior. 

PACS numbers: 05.45.-a,68.35.Rh,89.65.-s,89.75.Da 
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I. INTRODUCTION 

In the past few years there has been an increasing in- 
terest among theoretical physicists in complex phenom- 
ena occurring in fields that are far apart from the tra- 
ditional realm of Physics like Social and Economic sci- 
ences 0, 01 S 0- The main reason being that social 
and economic systems often exhibit many instances of 
complex dynamics, including self-organization, pattern 
formation, synchronization and phase transition-like phe- 
nomena that closely resemble those observed in nonequi- 
librium physical systems 0, S S ■ Physicists approach 
to these systems usually provides insights into the ba- 
sic ingredients that should be included in simple models 
in order to obtain the dynamics observed. Although it 
is clear that Physics inspired models of socioeconomic 
phenomena are often very simplistic views of very com- 
plicated systems, the aim is to show how complex macro- 
scopic dynamics might arise from rather simple rules op- 
erating at the 'microscopic' level of individual agents and 
their mutual interactions. 

In this paper we consider a very simple model of in- 
novation propagation dynamics in an economic system 
formed by agents ITol | . The aim is to describe in a sim- 
ple way the adoption of innovations that occurs among in- 
dustries, firms or individuals. Once a brand new product 
appears in the market, the agents should decide whether 
or not they will incorporate the new technology. Adopt- 
ing the new technology (in the form of a software, device, 
gadget, etc) has a cost, but at the same time it may im- 
prove business performance in the case of firms, or may 
level off life quality for individuals. Innovations are re- 
garded here in a broad sense and stand for any device 
or tool. For instance, a firm can decide to incorporate 
WWW technology by creating or revamping its WWW 
page, or going into e-commerce for the first time. A lay- 
man observation is that, if not always, in most cases, 
when the new technology actually improves performance 



its use will spread all over. 

In this approach two main mechanisms for the prop- 
agation of innovations are considered. Firstly, external 
pressure can push an agent to adopt an innovation. This 
mechanism intends to mimic exogenous influence, such 
as advertising, and is independent of the network struc- 
ture. Secondly, there is the interaction among agents, 
which depends on the underlying network structure and 
is introduced in the model by considering local coupling 
rules. A single tunable parameter C, which is fixed and 
the same for all agents, accounts for the agents' resis- 
tance to change and controls the dynamical behavior of 
the system. In earlier studies, some of us have already 
focused on the several outcomes of the model in the so- 
cial and economic context I9L ITTI IT2]. From the economic 
point of view, the main result is that the system presents 
an optimal behavior for an intermediate value of C, and 
that this can be quantified with a macroscopic observ- 
able. This feature is closely related with the statisti- 
cal properties of the profile of technological levels of the 
agents and its dynamical evolution. A proper character- 
ization of these properties can be done with the tools of 
statistical mechanics and it is the main aim of this work. 

In this paper we show that this model can be inter- 
preted as a surface growth model. Such interpretation 
allows us to analyze the dynamical behavior of the model 
as a kinetic roughening process akin to other nonequilib- 
rium surface growth systems. We find that the model 
exhibits a continuous phase transition between a rough 
and a flat phase at a critical value Cth of the control 
parameter. We focus on the scaling properties at the 
threshold in order to determine the critical exponents at 
the transition. By defining a convenient order param- 
eter and studying its finite-size scaling properties near 
criticality we are able to show that the horizontal corre- 
lation length diverges as ^ ~ \C — Cth\~'' , where v « 2.5. 
Close to the threshold, relaxation dynamics to the sta- 
tionary regime is characterized by diverging correlation 
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times T ~ where z « 0.57 is the dynamic exponent. 
The existence of a nonequihbrium roughening transition 
in a 1+1 dimension model makes it interesting also for 
statistical mechanics. It is known that phase transitions 
in nonequihbrium 1 + 1 dimensional systems are usually 
associated with systems with absorbing states In 
this case, the number of absorbing states and symmetries 
among them determine the universality class to which a 
particular system belongs to. Thus, it is of great interest 
to find models far from equilibrium which do not possess 
absorbing states but still display a phase transition. As 
we will see below, our model lacks absorbing states and 
the measured critical exponents suggest that this model 
belongs to a new universality class. Finally, we discuss 
the physical mechanisms behind the critical transition in 
this model. 
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FIG. 1: Snapshots of the profile of a system with = 1024 
for C = 0.5, C = 1.0, C = 2.0, and C = 5.0. 



II. THE MODEL 

We consider iV agents placed at the sites of a one- 
dimensional lattice with periodic boundary conditions. 
Each site (agent) i is characterized by a real variable hi. 
In general, we can consider this quantity as a character- 
istic of a given individual that other agents might want 
to imitate. When an agent has adopted a new feature 
(innovation), her neighbors become aware of the change 
and balance their interest (quantified as hi — hj) with 
their resistance to change C to decide if they would like 
to imitate this change. In this way C controls the mech- 
anism of imitation. This parameter is constant and the 
same for all the agents in the system. 

The system is updated as follows [lll |: 

1. At each time step an agent hi is randomly selected 
and 

hi^h^^ A, (1) 



present snapshots of the surface profile in the stationary 
regime for four different values of C. 

The time scale separation -namely, slow driving versus 
fast relaxation in the form of avalanches of activity- is 
similar to that occurring in self-organized critical (SOC) 
systems and dynamically drives the system towards a 
stationary state jl^]. We will see below that, at variance 
with most SOC systems, two different stable phases are 
possible: an ordered (flat) phase and a disordered (rough) 
phase with scale invariant properties. For small C, the 
driving process easily triggers avalanches that cover the 
whole system, leading to a uniform advance and a flat 
phase. On the contrary, for large C, there are almost no 
avalanches, and the system advances mostly due to the 
random updates, thus presenting an extremely heteroge- 
neous and rough profile. For intermediate values of C one 
can clearly see the presence of large avalanches and new 
updates. In fact, in the intermediate regime one can find 
the optimal growth regime in which the agents reach a 
given average level with a minimum number of upgrades 



where A is a random variable uniformly distributed 
in [0,1] [2^. The driving process accounts for the 
external pressure that may lead an individual to 
spontaneously update by adopting a new technol- 
ogy. This mechanism keeps the system out of equi- 
librium. 

2. The agents j S r(i), being the set of nearest- 
neighbours of agent z, upgrade if hi — hj > C. If the 
latter is satisfied, agent j imitates agent i by setting 
hj — hi . In this way the information of an update 
may spread beyond the neighbors of the originally 
perturbed site. This procedure is repeated until no 
one else wants to change, concluding an avalanche 
of imitation events. We thus assume that the time 
scale of the imitation process is much shorter than 
the one corresponding to the external driving. 

Starting from a flat initial condition, hi — for all i, 
the system evolves to a stationary regime. In Fig. ^we 



III. CRITICAL ROUGHENING TRANSITION 

A. Stationary regime 

In order to characterize the different morphological 
phases we have performed extensive numerical simula- 
tions of the model. The fluctuations of the profile height 
are measured by means of the global interface width (15| , 
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Where () stands for average over noise realizations. At 
each time step the mean height value 
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FIG. 2: Time evolution of the width W{t) for four different 
system sizes (iV = 512, N = 1024, N = 2048, and iV = 4096) 
when C — 0.5 and C — 2.0. Results correspond to averages 
over 500 realizations of the noise. The inset shows that Waat 
scales with system size. 



is also calculated. It is important to stress here that 
time is always measured in the external driving temporal 
scale, so that one time step t corresponds to an external 
update. As a consequence, the number of agents that 
change their state may vary from a single one (which 
changes from h to h + A) to any number of agents in the 
system if the update generates an avalanche. 

In the following we report on the behavior of the width 
in the two different phases. On the top panel of Fig. [21 
we show the behavior of W{N,t) for C — 0.5. The sat- 
uration value does not depend on the system size, which 
indicates that the system is in the smooth phase. On 
the bottom panel of Fig. |21 we show the numerical re- 
sults in the rough phase, for C = 2.0. In this case the 
saturation value Wsat (N) scales with the system size, as 
is shown in the inset. We find that in the rough phase 
the height fluctuations seem to fit reasonable well with 
a scaling as Wsat{N) ^ N'^-^^, which actually cannot 
be distinguished from a possibly logarithmic dependence. 
These results strongly suggest the presence of a roughen- 
ing transition. 

In order to study the critical behavior at the transi- 
tion threshold we introduce a convenient order parame- 
ter. When looking at the profile snapshots in the station- 
ary regime in Fig. ^ one can easily notice the presence 
of large plateaux, i.e. finite connected regions of agents 
that have the same height. The size of these flat regions 
decreases as C grows, since for C — > cxa the model has to 
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FIG. 3: Order parameter Mstat vs. C for four different system 
sizes, iV = 256, = 512, N = 1024, and iV = 2048. The 
points correspond to an average over 250 realizations of the 
noise. 



become equivalent to the random deposition model [15j. 
We found that the size of the largest plateau can be used 
as an order parameter. In the following, we shall call M 
to the size of the largest plateau in the system, normal- 
ized by the system size N. In this way, a completely flat 
profile corresponds to M — 1. We have also tried other 
common choices, as the often used Ising-like magnetiza- 
tion (1/N) and its variations How- 
ever, we found that our election has better scaling prop- 
erties for this particular case, since it takes into account 
the singular behavior of the flat phase in this model. 

Starting from a flat initial condition M{t) evolves un- 
til it reaches a stationary value. In Fig. Owe show the 
behavior of the stationary value of the order parameter 
Mstat{N,C) vs. C for four different system sizes. The 
order parameter allows us to distinguish the two phases 
discussed above. Note that for small values of the con- 
trol parameter, the system gets ordered, implying a flat 
phase. On the contrary, the stationary value of the order 
parameter goes to zero for large values of C as the system 
size becomes larger. Critical behavior is expected close 
to the threshold Cth and, as usual, it can be studied nu- 
merically by finite-size scaling techniques [l^ IT^ [T§| as 
follows. For any value of the control parameter C, there 
exists a horizontal correlation length ^, which diverges as 
f ~ e~'^ when the distance to the critical threshold goes 
to zero e = \C — Cth\ ^ 0. In finite systems this actually 
occurs for values of C close to, but not exactly at, the 
threshold since the finite-size critical behavior is encoun- 
tered as long as ^ '--^ TV, or equivalently when e ~ N~^/'^ . 
Close to the threshold, e ^ 0, for sufficiently large val- 
ues of the system size, Mgtat converges to a finite value 
obeying Mstat{N, e) ^ . Just at the critical point e = 
we expect the order parameter to decay as a power-law 
with the system size 



MstatiN,e = 0)^N 



(4) 



In Fig. 21 we plot our numerical results for Mstat {N, e) 
vs. N for different values of the distance to the threshold 
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FIG. 4: Order parameter Mstat vs. for four different values 
of the parameter C = 0.90, C = 0.95, C = 1.00 and C = 1.05. 
A power law decay Mstat ~ A'^"" *'' is observed for C — 1.0. 
Results correspond to averages over 500 realizations. 
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FIG. 5: Data collapse of the order parameter Mstat as given in 
Eq. Results correspond to averages over 500 realizations. 



e. Only for C — Cth a power-law with the system size 
can be obtained and the slope of the straight line in a 
log-log plot gives an estimation of the ratio P/v = 0.44 ± 
0.05 between critical exponents. We can thus identify 
Cth = 1-0 ± 0.1 with the critical point. After having 
determined the critical point, numerical data for different 
system sizes can be cast in the finite-size scaling ansatz 



(5) 



where the scaling function g{y) ^ const for y ^ 1, and 
9[y) ~ 2/^ if y ^ 1. In Fig. Owe plot a data collapse 
that allows us to determine the values of the exponents 
l/v = 0.40 ±0.05 and 13/v = 0.44 ±0.05. From these, we 
then have /3 ~ 1.10 and v ~ 2.50. 



B. Dynamics 

Since the model is out of equilibrium our study is com- 
pleted next with an analysis of the dynamic behavior, 
which contains much information about the universality 
of the roughening transition. In Fig. El we show the 



FIG. 6: Order parameter dynamics for three values of C: 
C = 0.90 in the smooth phase, the critical value C = 1.00 and 
C = 2.00 in the rough phase. Two system sizes, A'^ = 2048 and 
N = 4096, are represented to better appreciate the deviations 
from power-law behavior, indicated with a dashed line, for the 
values outside the critical region. The curves correspond to 
an average over 500 realizations. 



temporal behavior of the order parameter M{t^ N, C) for 
three different values of C (above, below and at the crit- 
ical threshold). Again, only at the critical point we may 
expect to find a power-law decay M{t, N,e = 0) ^ i-P/^t ^ 
where vt is the exponent associated with the diverging 
correlation time r ~ e""^' as e — > 0. The correlation time 
corresponds to the typical time that correlations survive 
in the system and is given by r ^ where z is the dy- 
namic exponent. The three exponents are related by the 
usual scaling relation z = vt/v provided dynamic scaling 
holds. In Fig. Elwe can see that only at the critical point 
a power-law behavior is observed, while deviations occur 
for C ^ Cth- The fit to a straight line in a log-log plot, 
as shown Fig. [7| leads to a determination of the ratio 
P/vt = 0.77 ± 0.05. One can write the dynamic scaling 
ansatz 



(6) 



for the order parameter, which at the critical point, e = 0, 
reads 



M{t, N,e = 0) = N-^'^f [t/N'') 



(7) 



where the scaling function f{u) ^ const for u ^ \ and 
f{u) - u-P/''' for M < 1. We can then use the values 
of the exponents just obtained to collapse our data as 
shown in Fig. [Slwith exponents vt ~ 1.43 and z ~ 0.57. 

Nonequilibrium phase transitions have been mostly 
related to the universality class of directed percolation 
(DP), with very few exceptions p^ . In particular, there 
are many examples of roughening transitions far from 
equilibrium that have been linked to DP, examples in- 
clude polynuclear growth models , solid-on-solid mod- 
els with evaporation at the edges of terraces 0| and the 
fungal growth model ^-^^ these systems, the DP 

process emerges at a particular reference height of the in- 
terface. In this case, the critical exponents characterizing 
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FIG. 7: Order parameter dynamics for six increasing values 
of the system size A'^, from top to bottom, 2® to 2^*. A power- 
law behavior M ~ t~°'''''^ is observed in the transient regime. 
Results correspond to averages over 500 realizations. 



FIG. 9: Fraction of sites / that can generate an avalanche 
with a single update as a function of C for three different 
system sizes, N = 100, N = 200 and N = 400. Curves are 
averaged over 1000 realizations. 
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FIG. 8: Dynamical data collapse of the order parameter at 
the critical point, as given by Eq. Q. The exponents used 
correspond to the ones obtained above, fi/v ^ 0.44 and z ~ 
0.57. 



the roughening transition can be obtained from those of 
DP, which in 1+1 dimension are given hy v — — 1.10 

= 1.73 



for the correlation length exponent, vt = 

for the time correlation exponent, and z = z^^ — 1.58 
the dynamic exponent Our results clearly suggest 

that the roughening transition occurring in the innova- 
tion propagation model does not belong to the DP class. 

The relation of many nonequilibrium critical models 
to DP has led to the proposal of the conjecture due to 
Janssen and Grassberger ,10- .2Q |. which states that a 
model belongs to DP under the following assumptions 
El: 

1. The model displays a continuous phase transition 
from a fluctuating active phase into a unique ab- 
sorbing state. 

2. The transition is characterized by a positive one- 
component order parameter. 

3. The dynamic rules involve only short-range inter- 
actions. 



4. Finally, the system has no special attributes like 
additional symmetries or quenched randomness. 

Any model satisfying all above four conditions has been 
found to belong to DP universality class, with no excep- 
tion to date. However, it is known that at least some of 
the above DP conditions can be relaxed. In fact, there 
are a few examples of systems that, despite exhibiting 
no absorb ing states [la, or having quenched dis- 

order I21L [2^ also display nonequilibrium phase transi- 
tions that belong to the DP universality class. Our model 
does not have absorbing states, since in both the rough 
and the flat phase the interface keeps fluctuating. Also 
and perhaps most importantly, interaction is not short- 
ranged, because of the avalanches of activity that give 
rise to nonlocal effects with finite probability. Their in- 
fluence in the dynamics is reflected by the extremly low 
value of the dynamic exponent, z = 0.57 < 2, signature 
of a highly super-diffusive behavior. It appears that this 
nonlocal interaction mechanism is the responsible for the 
deviation of the DP critical behavior. 

We believe that the transition takes place exactly at 
C — 1.0. This is directly related to the dynamical evo- 
lution rules of the model. We have deflned the external 
driving by choosing a random number from a uniform 
distribution in [0,1]. As a consequence, for C < 1.0, a 
random update on any site can generate an avalanche. 
On the other hand, for C > 1.0, only a small fraction of 
sites will be able to generate an avalanche with a single 
update. In order to quantify this effect we have studied 
the fraction of sites which can generate an avalanche with 
a single update. A site i with this property will satisfy 



h,±i + 1>C 



(8) 



In Fig. Owe present the fraction of sites / which are able 
to generate an avalanche as a function of C. The figure 
clearly shows that this fraction remains close to one for 
C < 1.0 and drops abruptly to a small value for C > 1.0. 
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IV. CONCLUSIONS 

In this work we have studied a simple model of inno- 
vation propagation dynamics in an economic system as 
a surface growth model. This has allowed us to char- 
acterize different morphological phases and also to an- 
alyze the dynamical behavior of the model as a kinetic 
roughening process. We have characterized a roughening 
transition and determined its critical exponents by finite- 
size scaling techniques. The value of the exponents do 
not coincide with known universality classes. We believe 
that the avalanche driven dynamics with its long-range 
effects is the reason why this model does not belong to 
the DP universality class. We have also presented a pos- 
sible mechanism for the transition occurring exactly at 

CtK = 1. 
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